Collisional decay of ®^Rb Feshbach molecules at 1005.8 G 
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We present measurements of the loss-rate coefficients Kam and K^m caused by inelastic atom- 
molecule and molecule-molecule collisions. A thermal cloud of atomic ^'^Rb is prepared in an optical 
dipole trap. A magnetic field is ramped across the Feshbach resonance at 1007.4 G. This associates 
atom pairs to molecules. A measurement of the molecule loss at 1005.8 G yields Kam = 2 x 
10~^° era? /s. Additionally, the atoms can be removed with blast light. In this case, the measured 



molecule loss yields K„ 



■ 3 X 10""' cm7s. 



Diatomic molecules associated from ultracold atomic 
gases using Feshbach resonances are in a highly-excited 
ro-vibrational state. An inelastic collision with another 
atom or molecule can lead to a vibrational de-excitation 
of the molecule. The difference in binding energy is re- 
leased as kinetic energy in the relative motion of the 
molecule and the collision partner. This energy is typi- 
cally much larger than the trap depth, so that both par- 
ticles escape from the trap. The first experiments on 
association of molecules with Feshbach resonances [1-8] 
showed that the lifetimes of molecules made from bosonic 
atoms differ largely from the lifetimes of molecules made 
from fermionic atoms. Measurements in bosonic systems, 
^^Na and ^^'^Cs, revealed loss-rate coefficients of typically 
5 X 10"^^ cm'^/s [9, 10]. In fermionic systems, ^Li and 
^"^K, the loss-rate coefficients far away from the Feshbach 
resonance are similar, but they can be suppressed by or- 
ders of magnitude when holding the magnetic field close 
to the Feshbach resonance [5, 11]. An explanation based 
on the Pauli exclusion principle for fermions was put for- 
ward [12, 13]. Recent experiments showed that the short 
lifetimes in the bosonic species can be circumvented in 
an optical lattice [14, 15]. 

Here we present a measurement of the loss-rate coef- 
ficients Kam and Kmm toT iuclastic atom-molecule and 
molecule-molecule collisions, respectively. The molecules 
are associated from atomic ®^Rb using the Feshbach res- 
onance at 1007.4 G [16] with a width of AB = 0.2 G 
[17, 18]. The loss measurements are performed at 
1005.8 G. In one measurement, atoms and molecules are 
in the trap simultaneously. The loss in this measurement 
is dominated by inelastic atom-molecule collisions and 
reveals Kam- In another measurement, remaining atoms 
are removed from the trap using blast light [7, 14] af- 
ter associating the molecules. The loss in this measure- 
ment is dominated by inelastic molecule-molecule colli- 
sions and reveals Kmm- Both measurements are per- 
formed in thermal clouds. The value of Kmm for a quan- 
tum degenerate cloud would be half as large [19]. 

In our measurements, the molecules are associated 
from atoms in their absolute ground state. Spontaneous 
dissociation of Feshbach molecules into unbound atom 
pairs with lower-lying spin states as observed in *^Rb 
is therefore impossible [20]. A previous measurement in 



® ''Rb using photo-associated molecules set an upper limit 
of Kam < 8 X 10~^^ cm^/s [21]. This limit is not ap- 
plicable in the present experiment, because a different 
ro-vibrational state is investigated and the experiment is 
performed at a very different magnetic field. 

A recent theoretical model for the state used here pre- 
dicts Kam = 3 X lO-i"^ cmVs at 1005.8 G [22]. Previous 
models [23, 24] are only applicable if the magnetic field is 
less than AB away from the Feshbach resonance, which 
is not the case in the present experiment. 

The experimental setup [16] is designed for the prepa- 
ration of a BEG of ^^Rb atoms in a magnetic trap. 
For the present experiment, the radio- frequency induced 
evaporation is stopped near the critical temperature Tc 
of the phase transition to BEG. The atoms are then 
transferred into a crossed-beam optical dipole trap as de- 
scribed in Ref. [18]. The dipole trap is operated at some- 
what higher laser power than in Ref. [18]. The measured 
trap frequencies are {uJx,i^y, i^z) = 2nx (95, 154, 200) Hz. 

After transfer into the dipole trap, a magnetic field 
oi B = 1007.6 G is turned on rapidly and the spin 
of the atoms is transferred [17] to the hyperfine state 
[F = 1, mF = 1). Next, B is ramped across the Feshbach 
resonance at a rate of 0.4 G/ms. This associates atom 
pairs to molecules as described in Refs. [3, 25]. As soon as 
molecules are forming during the ramp, they can undergo 
inelastic collisions. For molecules made from bosonic 
atoms, the inelastic collision rates are enhanced near the 
Feshbach resonance [7] . Hence, the molecule number can 
be maximized by jumping the magnetic field away from 
the resonance as fast as possible, once the molecules are 
created [26]. To this end, B is jumped from 1007.35 G 
to 1005.8 G. We find experimentally that this combina- 
tion of ramp speed and start point for the magnetic-field 
jump produces the maximum molecule number. 

Immediately after jumping the field to 1005.8 G, re- 
maining atoms can be removed from the trap by applying 
blast light for 0.3 ms as described in the appendix. This 
is followed by a variable hold time in the trap. This time 
is scanned in the loss measurements described below. At 
the end of this hold time, the trap is switched off. Right 
after release from the trap, the molecules are separated 
from remaining atoms using the Stern-Gerlach effect by 
applying a magnetic-field gradient of 120 G/cm for 1 ms. 



2 



Immediately after this, the magnetic field is jumped to 
1006.9 G and subsequently the molecules are dissociated 
into unbound atom pairs by ramping the magnetic field 
back across the Feshbach resonance to 1007.7 G at a rate 
of 0.8 G/ms. At the end of this ramp, the magnetic field 
is switched off rapidly. Finally, 7 ms after release from 
the trap, an absorption image is taken. 

The molecule number decays as a function of hold time 
between association of the molecules and release from 
the trap. The loss of molecules from the trap can be 
described by the rate equation 
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where tIq and the particle densities of atoms and 

molecules and Kam and Kmm are the loss-rate coeffi- 
cients caused by inelastic atom-molecule and molecule- 
molecule collisions, respectively. K^^ represents molecule 
loss mechanisms which do not rely on collisions with 
other cold atoms or molecules. Such loss could be caused 
by background gas collisions, photo-dissociation by the 
dipole-trap light, or spontaneous decay into lower ro- 
vibrational levels. Our experimental results show that 
Km is negligible. 

The loss of atoms during the hold time is also found to 
be negligible. This is because the atom number is either 
zero or much higher than the molecule number, so that 
inelastic atom-moleculc collisions can only lead to loss of 
a small fraction of the atoms. 

Volume integration of Eq. (1) yields 
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where Na and Nm are the total number of atoms and 
molecules, respectively, and where we abbreviated 
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N,N„ 



Ij.Tlm. d X 



(3) 



for i equal to a or m for atoms or molecules, respectively. 
Vim is an effective volume and depends on the shape of 

the cloud, but not on the particle number. Note that 
{n-m) = ^m/ymm is often referred to as the average den- 
sity. 

Assuming that Na, Vam, and Vmm are time indepen- 
dent, the rate equation (2) can be integrated analytically 
using standard methods, yielding 



Nm{t) 
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(4) 



where /3 = Kmm/Vmm and T = Km + KamNa/Vam and 
A^o is the molecule number at t = 0. 

Fig. 1 shows experimental results for the molecule loss 
obtained after blasting away the atoms. In this measure- 
ment, r = K„i because Na = 0. An unconstrained fit 
of Eq. (4) to the data yields a slightly negative value for 
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FIG. 1: Loss of molecules after blasting away the atoms. The 
solid line shows a fit of Eq. (4) to the data, yielding Kmm- 
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FIG. 2: Loss of molecules in the presence of the atoms. The 
solid line shows a fit of Eq. (4) to the data, yielding Kam- 
Clearly, the loss here is much fcister than in Fig. 1. 



r which is unphysical. We therefore fix F = and ob- 
tain /J = (8 ± 1) X 10~^/s from the fit. The error bar is 
statistical. 

Fig. 2 shows experimental results without blasting 
away the atoms. Obviously, the presence of the atoms 
reduces the molecule lifetime substantially. An uncon- 
strained fit of Eq. (4) to the data yields a slightly nega- 
tive value for (3 which is unphysical. We therefore fix /3 
to the value from Fig. 1 and obtain F = (0.9 ± 0.1)/ms 
from the fit. Again, the error bar is statistical. 

In order to extract the loss-rate coefficients Kmm and 
Kam, the effective vohimcs Vmm and Vam must be deter- 
mined. This requires knowledge of the density distribu- 
tions of the atomic and molecular cloud. This is a delicate 
issue, because the small cloud size makes direct measure- 
ments of the spatial distributions in the trap very hard. 
Theoretical modeling is also hard. Even the distributions 
at the beginning of the loss measurements are difficult to 
model, because the time evolution during molecule asso- 
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ciation is nontrivial. This is because, first, atomic pair 
correlations are crucial for a realistic treatment of the as- 
sociation process [27]. Second, the loss rate coefficients 
a-nd Kam become relevant as soon as molecules 
start to form during the ramp. Third, these coefficients 
exhibit an unknown but probably strong magnetic-field 
dependence near the Fcshbach resonance. 

During the loss measurement, anharmonicities in the 
trap potential and elastic collisions between particles 
tend to randomize the motion leading towards a thermal 
distribution. Inelastic collisions, however, remove parti- 
cles preferentially at the trap center, where the density 
is highest. The elastic scattering cross sections for the 
molecules are unknown. Hence, it is not clear which pro- 
cess dominates at what stage of the loss measurement. 
The evolution of the density distributions during the loss 
measurements is therefore a complex process with un- 
known parameters and unknown initial conditions. 

In order to obtain an estimate for the effective vol- 
umes Vam and Vmm ) we assume that particles of the same 
species (atoms or molecules) are in thermal equilibrium. 
Our model does allow for a temperature difference be- 
tween atoms and molecules. With this approximation, 
the spatial density distributions Ua and Gaussian 
and Eq. (3) yields the effective volume for species i 

Vim = {2n f' II ^al,+al^,, (5) 
fe=i 

where the index k refers to the three directions in space 
and the one-dimensional (ID) root-mean-square (rms) 

radii of the Gaussians are ai^k = ^kBTi/(m.iUjf ^,) . The 

mass and temperature of species i are labeled rrii and 
Ti, respectively. The polarizability and the mass of a 
molecule arc twice as large as for one atom, so that the 
trap frequencies for atoms and molecules are identical, 

i.e. Ua^k = <JJm,k- 

The temperature of the atoms and molecules is deter- 
mined from time-of-fiight measurements. For the atoms, 
the cloud size evolves as 

OaAt) = ^J<^lm + <at\ (6) 

where the trap is switched off at t = and cfy^i = 
^JkBTi/rrii is the ID rms-velocity of species i, which is 
independent of the spatial direction if the species is in 
thermal equilibrium. For the molecules, extra kinetic en- 
ergy is added in the dissociation process [9, 18], so that 

= \[^l^M + <m^+<A. ■ (7) 

Again, i = is chosen at the time of release from the 
trap. f7v,dis refiects the extra kinetic energy released in 
the dissociation, t-^cm is the remaining time of flight be- 
tween dissociation and detection. In order to determine 



Tjn oc CT^ ^ , we scan the time t between release and de- 
tection in such a way remains fixed. A fit of 
Eq. (7) is then equivalent to a fit of Eq. (6) with a mod- 
ified value for <t(0). Hence, the extracted temperature is 
insensitive to the dissociation heating. 

The atomic cloud before molecule association typically 
contains 3.6 x 10^ atoms at a temperature of 0.5 /iK very 
close to Tc- The cloud is almost purely thermal with only 
6 X 10^ BEG atoms. The molecule association is accom- 
panied by noticeable heating and substantial loss. The 
remaining atomic cloud contains Na = 1.9 x 10^ atoms 
at a temperature of = 1.0 fiK. There is no BEG in the 
remaining atomic cloud. The molecular cloud has a tem- 
perature of Tm = 1.5 fiK resulting in a peak phase-space 
density of ^ 10~^ for the molecules. The center-of-mass 
motion of an atom pair is unchanged in the association 
and in a dilute thermal cloud the atomic pair correlations 
are uncorrelated from the center-of-mass motion of the 
pairs. Hence, our experiment should produce molecules 
with the same temperature as the initial atoms. This 
agrees reasonably with our measurements. The mea- 
sured values of Na, Ta, and Tm vary by less than 10% 
during the loss measurement. This justifies the assump- 
tion that they are time- independent, which was used to 
derive Eq.(4). 

The above values for and yield effective volumes 
of Vmm = 3.8 X 10"^ cm^ and Km = 4.6 x 10"^ cm^. 
The resulting loss-rate coefficients are Kmm = 3 x 
10"^° cm^/s and Kam = 2 x 10"^° cm^/s. Statistical 
errors on the rate coefficients are negligible compared to 
systematic errors. 

The dominant systematic error arises from the prob- 
lematic assumption that the clouds are in thermal equi- 
librium. For reasons discussed above, it is hard to mea- 
sure or model the temporal evolution of the density distri- 
butions. On one hand, the association process preferen- 
tially populates the trap center, because the association 
is less efficient at low atomic density. On the other hand, 
inelastic collisions preferentially deplete the trap center. 
Both effects can lead to a misestimation of Vam but the 
trends go into opposite directions. It is not clear which 
of the effects dominates. The resulting systematic error 
is hard to quantify. We speculate that a factor of 3 seems 
possible for Kam as well as Kmm- 

In conclusion, we measured the loss-rate coefficients 
Kam and K^m for Feshbach molecules in ^''Rb at 
1005.8 G. These results yield valuable input for theo- 
retical models of the three-atom and four-atom system. 
The measured value for Kam agrees reasonably with the 
prediction Kam = 3 x 10"!° cm^/s of Ref. [22]. 

Appendix: Blast light 

After associating molecules, remaining atoms can be 
pushed out of the trap using the radiation pressure of 
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applied laser light [7, 14]. The dipole trap used here is 
much deeper than the photon recoil energy, so that the 
blast light must drive a closed cycling transition, which is 
not possible when starting from the state \F = l.mp = 

1) . Hence, a first light field (the "pump light") optically 
pumps the atoms from |F = 1, mp = 1) to |F = 2, mp = 

2) . A second light field (the "cycling light") then drives 
a cycling transition. We only drive transitions on the 
5^<S'i/2 5^P3/2 resonance line at 780 nm. 

As the blast light is applied while the atoms are in 
a magnetic field of 1005.8 G, some effort is needed to 
obtain light stabilized at the required frequencies. Two 
frequency-stabilized lasers arc already in use in the exper- 
iment for operating the magneto-optical trap (MOT): one 
is locked close to the MOT transition \F = 2) ^ \F' = 3) 
and another laser is locked to the rcpump transition 
\F = 1) <-> \F' = 2). Quantum numbers with a prime 
refer to excited states. 

The cycling laser resonantly drives the (t"*" transition 
|F = 2,mF = 2) ^ \m'j = 3/2, = 3) at 1005.8 G, 
which is 1405 MHz blue detuned from the B = 
MOT transition. Note that at 1005.8 G, the relevant 
ground states are characterized by good quantum num- 
bers F,mF, while the excited states are characterized by 
good quantum numbers m'j,m'p. A beat lock stabilizes 
the frequency of the cycling laser relative to the MOT 
laser. 

In order to deplete the state \F = l,mp = 1), we 
use pump light that is 62 MHz red detuned from the 
B = repump transition. The light is obtained from the 
rcpump laser using an acousto-optical modulator (AOM) 
and resonantly drives the tt transition |F = l,mp = 
1) ^ \m'j = 3/2, m'p = 1) at 1005.8 G. 

From the excited state populated by this pump light, 
atoms can decay back to the initial state, into the de- 
sired state \F = 2,mp = 2) or into the undesired state 
\F = 2,nip = 1). The experiment shows, that the 
branching ratio for decay into the undesired state is only 
a few percent. Still, we use a third light field to de- 
plete this state. Wc use another AOM to obtain light 
that is 87 MHz red detuned from the B = MOT tran- 
sition. This light resonantly drives the cr+ transition 
\F = 2,mp = l) ^ \m'j = 3/2, m'p = 2) at 1005.8 G. 

Note that the excited state \m'j = 3/2, m'p = 1) can- 
not decay into ground states with mp = Q, because at 
1005.8 G the quantum number mj = 3/2 is conserved 
during the decay and the J = 1/2 ground state does not 
have sub-states with mj = —3/2. 

All three blast fields are operated a factor of ten or 
more above saturation intensity. They are on simultane- 
ously. The trap depth is a few microkelvin, so that ap- 
proximately five directed photon recoil momenta should 
add sufficient kinetic energy for an atom to leave the trap. 
This momentum should be accumulated after 0.2 /zs. The 
estimated acceleration is 10^ m/s^, so that the cloud ra- 
dius of 10 fim is estimated to be traversed in 15 /xs. But 



the experiment shows that the blast light needs to be on 
for 300 /zs, in order to remove all atoms from the trap. 
The origin of this discrepancy is unclear. 

Application of the three blast beams reduces the 
molecule number by approximately 30%. A recent exper- 
iment, performed after the measurements reported here, 
demonstrated that atoms can be removed by applying a 
microwave field and only the cycling laser [14]. This has 
the advantage that the loss of molecules during applica- 
tion of the blast light becomes negligible [15]. 
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